Angular dependence of magnetic properties in Ni nanowire arrays 
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The angular dependence of the remanence and coercivity of Ni nanowire arrays produced inside 
the pores of anodic alumina membranes has been studied. By comparing our analytical calculations 
with our measurements we conclude that the magnetization reversal in this array is driven by means 
of the nucleation and propagation of a transverse wall. A simple model based on an adapted Stoner- 
Wohlfarth model is used to explain the angular dependence of the coercivity. 

PACS numbers: 75.75.+a,75.10.-b,75.60. Jk 
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I. INTRODUCTION 

Arrays of magnetic nanowires have attracted consider- 
able interest due to their promising technological appli- 
cations [H-Q mainly in high-density information storage. 
Electrodeposition in polycarbonate and alumina mem- 
branes is the preferred technique used in the fabrication 
of these systems [EEl, leading to the production of highly 
ordered arrays of magnetic nanowires inside the pores 
of membranes ill. Different groups have investigated 
the role of magnetostatic interactions and the influence of 
the size of nanowires on the magnetic properties of these 
systems [ITl - fl4| . Recently Escrig et al [Tl| investigated 
the dependence of the coercivity of Ni nanowire arrays 
when the external magnetic field is applied parallel to 
the nanowire axis. However, for applications in devices 
that use longitudinal/perpendicular magnetic recording, 
it is important to investigate the angular dependence of 
the magnetization. It is expected that the direction of 
the external magnetic field with respect to the wire axis 
strongly influences the magnetic properties of the array. 

The properties of virtually all magnetic materials 
are controlled by domains separated by domain walls. 
Measurements on elongated magnetic nanostructures [l5j 
highlighted the importance of nucleation and propaga- 
tion of a domain wall between opposing magnetic do- 
mains in the magnetization reversal process. Three main 
idealized modes of magnetization reversal have been iden- 
tified and occur depending on the geometry of the par- 
ticles m m, mi. These mechanisms are known as co- 
herent rotation, C, with all the spins rotating simultane- 
ously; vortex wall, V, in which spins invert progressively 
via propagation of a vortex domain wall; and transverse 
wall, T, in which spins invert progressively via propaga- 
tion of a transverse domain wall. If a wire is thin enough, 
the exchange interaction forces the magnetization to be 
homogeneous through any radial cross section of the par- 
ticle (lcf . In this way a transverse wall is the preferred 
mode in thin ferromagnetic wires (diameters d < 60 nm) . 

To investigate the angular dependence of the rever- 
sal of the magnetization in nanoparticles, several models 



have been proposed. In particular, the Stoner-Wohlfarth 
model has been used to calculate the angular dependence 
of the coercivity when the reversal of the magnetization 
is driven by coherent rotation [19j . However, Landeros 
et al [13] have shown that coherent rotation is present 
only in very short particles, namely, when the length of 
the particles is similar to the domain wall width. On 
the other side, Aharoni [2(| calculated the angular de- 
pendence of the nucleation field in an ellipsoid when the 
reversal is driven by curling rotation. Escrig et al (2lj 
extended this expression for magnetic nanotubes. Fi- 
nally, Allende et al [22| presented an expression for the 
angular dependence of the coercivity in magnetic nan- 
otubes when the reversal of the magnetization is driven 
by a transverse domain wall. However, the angular de- 
pendence of the coercivity in magnetic nanowire arrays 
has not been modeled yet. Thus, usually experimental 
results for the angular dependence of the coercivity are 
interpreted using the Stoner-Wohlfarth model |19| . 

In this paper we present an analytical model that al- 
lows us to investigate the angular dependence of the coer- 
civity and remanence for Ni nanowire arrays, considering 
the different modes that can be present, as a function of 
wire geometry. Additionally, experimental data for the 
switching field of high-aspect ratio Ni nanowires will be 
compared with this analytical model. 



II. EXPERIMENTAL METHODS 

Our approach to the preparation of magnetic Ni 
nanowires arranged in hexagonally ordered, parallel ar- 
rays is based on the combination of two complementary 
aspects, namely, (i) the use of self-ordered anodic alu- 
minum oxide (AAO) as a porous template and (ii) the 
electrodeposition of Ni in the cylindrical pores. 

Anodic aluminum oxide is obtained from the electro- 
chemical oxidation of aluminum metal under high volt- 
age (usually 20-200 V) in aqueous acidic solutions [f|. 
Under certain proper sets of experimental conditions 
(acid nature and concentration, temperature, and ap- 
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plied voltage), the electrochemically generated layer of 
alumina displays a self-ordered porous structure. Cylin- 
drical pores of homogeneous diameter are thus obtained, 
with their long axis perpendicular to the plane of the 
alumina layer and ordered in a close-packed hexagonal 
arrangement. With our method [23j], first and second 
anodizations of Al have been carried out in 0.3M oxalic 
acid under 40 V at 20 °C. Subsequently, the pores were 
widened in a 0.085M H3PO4 solution at 37 °C, obtaining 
pores of d = 2R = 50 nm diameter and a center-to-center 
distance of D = 100 nm. 

Electrodeposition was performed at a constant poten- 
tial (dc electrodeposition at -1.0 V). To facilitate the elec- 
tric contact, a very thin Au-Pd layer was sputtered on one 
side of the membrane, followed by the electrodeposition 
[ill I24I ] of a thicker nickel layer to achieve full pore seal- 
ing. The potentiostatic condition allows us to have more 
precise control of the electrochemical reaction, and then 
more accurate control of the growth of the Ni nanowires 
111. 



The morphology and the single-crystal structure of the 
individual nanowires after dissolution of the template 
were subsequently investigated by scanning electron mi- 
croscopy (SEM) and transmission electron microscopy 
(TEM) using a JEOL 5900 LV, checking the high or- 
dering of the hexagonal arrays and the large aspect ra- 
tio. [ll| The chemical characterization of the nanowires 
was made by means of energy dispersive analysis of X- 
rays (ED AX), obtaining a clear peak corresponding to 
the emission spectrum of Ni. |llj The magnetic mea- 
surements were performed by a vibrating sample magne- 
tometer (VSM) with the applied field at different angles. 
The measurements were made at ambient temperature. 
Figure 1 illustrates the hysteresis curves for two samples 
with diameter d — 50 nm, lattice parameter D = 100 nm, 
and lengths (a) L — A /im and (b) 12 /jm, measured with 
the external field along and perpendicular to the axis of 
the wires. In this figure it is clearly seen that the easy 
axis of the nanowires is along the wire axis. 

In a series of Ni nanowire arrays of varying lengths (all 
other geometric parameters being kept constant) investi- 
gated by VSM magnetometry, we observed a significant 
dependence of the remanence (see Fig. 2) and coercivity 
(see Fig. 3) on the angle at which the external magnetic 
field is applied. In particular, the coercive field H c can be 
tuned between 100 and 720 Oe approximately by prop- 
erly adjusting 9. The 12 (im nanowires present a sharp 
decrease of the coercivity and remanence near 90° , while 
the 4 fim sample has a broader curve. 



III. ANALYTICAL MODEL 

In order to analyze our results for the angular depen- 
dence, we made calculations which lead us to obtain an- 
alytical expressions for both coercivity and remanence. 
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FIG. 1: Hysteresis curves at external field parallel and per- 
pendicular to the nanowire axis, for two arrays with (a) L = 
4 fim and (b) 12 /im, d — 50 nm, and D = 100 nm. 



A. Angular dependence of the remanence 

In an array the magnetization inside each nanowire 
is oriented parallel to the wire's axis (easy axis), due 
to the strong shape anisotropy. When the applied field 
is reduced to zero, at remanence, each wire presents its 
magnetization along the axis, but it is measured at an 
angle 9 with respect to the easy axis. When the field is 
reduced to zero, and due to the strong shape anisotropy, 
the remanence rapidly relaxes and points along the axis. 
Then, one can approximate the remanence of the wire 
by the projection of the magnetization on the direc- 
tion of measurement, i. e., Mr (9) = Mfl |cos6*|, with 
Mr = Mr (9 = 0) the remanence measured at 9 = 
0. The solid line in Fig. 2 illustrates this expression 
for Mr (9) /Mq together with our experimental results, 
showing a very good agreement, particularly for L = 12 
jj,m, for which the shape anisotropy is stronger. Then the 
angular dependence of the remanence can be described 
by Mr (9) = Mr\cos9\ for nanowires with a large as- 
pect ratio, i. e., with their easy axis along the wire axis. 
For nanowires in which the easy axis is perpendicular to 
the wire axis, the angular behavior of the remanence can 
be modeled as Mr \sin9\. Note that this expression de- 
fines only the angular behavior of the remanence, and 
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FIG. 2: Experimental results for the reduced remanence, 
Mr (6) /Mo (where Mo is the saturation magnetization) as a 
function of the angle 9 (dots), and analytical expression given 
by (Mr/Mo) |cost9|, using the experimental values of the in- 
trinsic remanence, Mr, obtained from the hysteresis loops of 
Fig. 1 (solid line). 



the reversal mechanisms, k = C, T, and V. The angular 
dependence of the nucleation for a coherent magnetiza- 
tion reversal was calculated by Stoner-Wohlfarth [19| and 
gives 
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where t — tans (0) and Mq is the saturation mag- 
netization. The demagnetizing factor of a wire along 
the z axis has been previously obtained [H, HI] and is 



given by N z (I) ■ 
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] is a hypergeometric function. 
For the T mode, the angular dependence of the co- 
ercivity can be studied by an adapted Stoner-Wohlfarth 
model [2l[|22| in which the width of the domain wall, wt, 
is used as the length of the region undergoing coherent 
rotation. Starting from the equations presented by Lan- 
deros et al [l7| we can calculate the width of the domain 
wall for the transverse mode as a function of the wire ge- 
ometry. The domain wall width increases with the wire's 
radius, and is about 50 nm larger than the wire's radii. 
Following this approach, 
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FIG. 3: Experimental data for the angular dependence of 
coercivity in Ni nanowire arrays, with length L — 4 fim and 
12 fim, d = 50 nm, and D — 100 nm. 
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It is important to mention that the expressions for (0) 
and (9) differ only by the length. In a coherent re- 
versal, L represents the total length of the wire, and then 
the coercivity varies with the length. However, in a trans- 
verse mode, and because u>t is almost independent of the 
length of the wire, as shown by Landeros et al [17| . the 
coercivity is also independent of the length. Then, when 
L > wt {L < wt), the T (C) mode will always exhibit a 
lower coercivity, independent of 8. 

As shown in the Stoner-Wohlfarth model [l!|, the nu- 
cleation field does not represent the coercivity in all cases. 
However, from the discussion presented on p. 21 in ref- 
erence [19], the coercivity for coherent and transverse 
reversal modes can be written as 
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< 9 < it/4 



2 Hn (T) (9 = 7r/4) - Hn (T) ir/4<9< ir/2 



does not determine the intrinsic remanence, Mr, which 
has to be obtained by considering the particular reversal 
mode, according to the geometry of the wire. 



B. Angular dependence of the coercivity 

In order to study the angular dependence of the re- 
versal modes we have developed analytical calculations 
that lead us to obtain the coercive field H 1 ? for each of 



In very short wires (L k, wt) the energy cost involved 
in the creation of a domain wall gives rise to a coherent 
reversal. 

The curling mode was proposed by Frei et al [27j and 
has been used to investigate magnetic switching in films 
28 1 and particles with different geometries, like spheres 
27j, prolate ellipsoids [13, H3], and cylinders How- 
ever, for simplicity, expressions for the nucleation field 
obtained using infinite cylinders are used in all cases. 
The angular dependence of the curling nucleation field in 
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a finite prolate spheroid was obtained by Aharoni[2C 
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The demagnetizing factor of a wire along the x axis 
has been obtained previously [25|, [26| and is given by 

N x (I) = \l 
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For a cylindrical geometry, 

Shtrikman et al [29| have obtained q 2 = I.O871". Here 
J p (z) and Y p (z) are Bessel functions of the first and 
second kind, respectively, and L x = y/2A/ n n M 2 is the 
exchange length. As pointed out by Aharoni [20|], for a 
prolate spheroid with 9 = 0, a jump of the magnetization 
at or near the curling nucleation field occurs. Therefore, 
the coercivity is quite close to the absolute value of the 
nucleation field. Then, we assumed here that — is a 
good approximation to the coercivity, , when the re- 
versal occurs in the V mode, as in other studies [29l. l3l|. 
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FIG. 4: Magnetization reversal modes in nanowires by coher- 
ent rotation mode and transverse reversal mode. 



IV. RESULTS AND DISCUSSION 

With the above expressions we can study the angular 
dependence of the coercivity for our Ni nanowire array. 
In our model the system reverses its magnetization by 
whichever mode opens an energetically accessible route 
first, that is, by the mode that exhibits the lowest coer- 
civity. However, for highly dynamic cases, the path of 
lowest coercivity might not be accessible due to preces- 
sion effects. By evaluating the coercivity for the different 
modes described in Section III we found the one which 
drives the reversal for each 9. Figure 5 illustrates our 
results for an isolated Ni nanowire with L = 12 fim. The 



dotted line represents the coercive field when the wire 
reverses its magnetization by a coherent rotation; the 
dashed line is the coercive field when the wire reverses 
its magnetization by the propagation of a vortex domain 
wall, and the solid line represents the coercive field when 
the wire reverses its magnetization by the propagation of 
a transverse domain wall. 
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FIG. 5: Figure 5: Coercivity for different angles between the 
applied magnetic field and the wire axis (from 0° to 90°). 

As depicted in Fig. 5, in a 50 nm diameter Ni nanowire 
the T mode is preferred for almost every 9 because it 
exhibits a lower coercivity. However, in the limit 9 ~ 90, 
the coercivity for the T and C modes are very close. 

Our results for Ni nanowires are combined in Fig. 
6. The solid line represents the coercivity obtained by 
means of analytical calculations for L = 12 /im. However, 
for almost every 9, differences with the results for L = 4 
fxm are less than 3 %, then at this scale, the line depicts 
results for both lengths. Squares and triangles represent 
the experimental values of the coercivity for L = 12 /im 
and L = 4 /zm, respectively. As shown in Fig. 5, the co- 
ercivity obtained from the Stoner-Wohlfarth model (co- 
herent rotation) is two to three times higher than the one 
obtained assuming a transverse domain wall, for almost 
all 9, and is much larger than the experimental results. 
Then, although the shape of the angular dependence of 
the coercivity is the same for coherent rotation and the 
transverse reversal mode, the values are significantly dif- 
ferent, showing that a transverse mode drives the reversal 
for each 9. 

However, the absolute values computed for the coer- 
civity are still different compared with the experimen- 
tal data. We ascribe such difference between calcula- 
tions and experimental results to the interaction of each 
wire with the stray fields produced by the array. When 
the distances between nanowires are smaller or compa- 
rable to their diameter, the magnetostatic interactions 
between nanowires are stronger, and consequently have 
an important influence on the magnetic properties of the 
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FIG. 6: Angular dependence of the coercive field measured 
experimentally and calculated analytically for Ni nanowires 
with 50 nm diameter. 



array. The interaction of each wire with the stray fields 
produced by the array strongly influences the coercivity 
[33 - l34| . Therefore, the small discrepancy between ex- 
periments and model can be regarded as the result of 
interactions within the wires in the array and size distri- 
bution, which are not included in our model. 

Finally, the results presented above may be general- 
ized. Figure 7 shows the angular dependence of the co- 
ercivity for an isolated Ni nanowire when it reverses its 
magnetization by the propagation of a transverse domain 
wall. We have considered different radii in order to in- 
vestigate the behavior of the geometry of the wire with 
the angular dependence of the coercivity. In the range of 
radii considered, we find that an increase in the radius R 
results in a decrease in the coercive field. 



suming that the relaxation process is very fast and then 
the magnetization orients rapidly along the easy axes. 
We have derived analytical expressions that allow us to 
obtain the coercivity when the wire reverses its magneti- 
zation by means of a coherent rotation, transverse rever- 
sal mode, and a vortex domain wall. For the wires stud- 
ied experimentally the magnetization is driven by means 
of the nucleation and propagation of a transverse wall. 
The stray field produced by the array, as a function of 
the angle of the external field, cannot be determined an- 
alytically. However, it is responsible for small differences 
between experiments and calculations. 
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FIG. 7: Angular dependence of the coercivity in an isolated 
Ni nanowire with different values of its radius. We have con- 
sidered a transverse reversal mode. 
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